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■ Abstract. Sharp upper and lower bounds are established for the modes of the Poisson 

I , distribution of order k. The lower bound established in this paper is better than the previously 

' established lower bound. In addition, for A; = 2, 3, 4, 5, a recent conjecture is presently proved 

I solving partially an open problem since 1983. 

^ . 

1. Introduction and Summary 

For any given positive integer k, denote by Nk the number of independent trials with 
i constant success probabihty p until the occurrence of the k^^ consecutive success, and set 

1 — p. For n > k, Philippou and Muwafi [13j derived the probability P{Nk = n) in terms 



00 



' of multinomial coefficients and noted that P{Nk = n\p = 1/2) = /^^^;j,_,_]^/2" where /n*^^ is the 

^ ■ n*^ Fibonacci number of order k. See also [3], [15], and [16]. Philippou et al. [12] showed that 

X^^j. P{Nk = n) = 1 and named the distribution of A^^ the geometric distribution of order 
k with parameter p, since for A; = 1 it reduces to the geometric distribution with parameter 
. p. Assuming that A''^ i, . . . , A''^ are independent random variables distributed as geometric of 

I order k with parameter p, and setting Y^^^ = ^k,j, the latter authors showed that 

^ P{Yk,r = y)=p'Y.(^'^"'^^'^'^^) i-r^-^'^ y = kr,kr + l,..., 

rn . ^ \ yi,- ■ ■ ,yk,r - 1 J p 

^ . where the summation is taken over all fc-tuples of non- negative integers yi,y2,- " iVk such 

that ui + 2y2 + • • • + ky^ = y — kr. They named the distribution of Y^^^ the negative binomial 
distribution of order k with parameters r and p, since for fc = 1 it reduces to the negative 
binomial distribution with the same parameters. Furthermore they showed that, if rq — )• 
A (A > 0) as r — )• oo and g — )• 0, then 



o 
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lim P{Yk,r - kr = x) = V e"^^ ^ '\ ^ = U^-A), x = 0,1,2,..., (1.1) 

where the summation is taken over all fc-tuples of non- negative integers xi,X2, • • • , x^ such 
that xi + 2x2 + • • • + kxk = x. They named the distribution with probability mass function 
fk{x; A) the Poisson distribution of order k with parameter A, since for /c = 1 it reduces to the 
Poisson distribution with parameter A. See also [I], |9], and [2]. 

Denote by nik^x the mode(s) of fk{x;X), i.e. the value(s) of x for which fk{x;X) attains 
its maximum. It is well known that rui^x = A or A — 1 if A G N, and mi^\ = [AJ if A N. 
Philippou [7] derived some properties of fk{x; A) and posed the problem of finding its mode(s) 
for k>2. See also [8] and [11] • 

Hirano et al. [5] presents several graphs of fk{x; A) for A € (0, 1) and 2 < A; < 8, and Luo 
[6] derived the following inequality 



mk,\>kXWkl ^— — k>l{X>0), (1.2) 
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which is sharp in the sense that nii^x = A — 1 for A G N. Recently, Phihppou and Saghafi |14] 
conjectured that, for A: > 2 and A G N, 

mk,x = Xk{k + l)/2-lk/2\, (1.3) 

where [u\ denotes the greatest integer not exceeding n G M. 

In this paper, we employ the probability generating function of the Poisson distribution 
of order k to improve the bound of Luo [6] and to also give an upper bound (see Theorem 
2.1). We then use Theorem 2.1 to prove the conjecture of Philippou and Saghafi [14] when 
k = 2,3, 4, 5, partially answering the open problem of Philippou [71 [HI [TT]. 

2. Main results 

In the present section, we state and prove the following two theorems. 

Theorem 2.1. For any integer k > 1 and real A > 0, the mode(s) of the Poisson distribution 
of order k satisfies (y) the inequalities 

[Xk{k + 1)/2J - M^±ll + 1 _ Sk,i < mfc,A < [Xk{k + 1)/2J , 
where 6k^i is the Kronecker delta. 

Theorem 2.2. For A G N and 2 < k < 5, the Poisson distribution of order k has a unique 
mode rrik^x = Xk{k + l)/2 - [A;/2j . 

For the proofs of the theorems we employ the probability generating function of the Poisson 
distribution of order k and some recurrences derived from it. We observe first that the left 
hand side inequality in Theorem 2.1 is sharp since, for A G N, mi^A = X — 1, the value of the 
lower bound for k = 1. The right hand side inequality is also sharp in the sense that there 
exist values of k and A for which ruf^ x = [Xk{k + l)/2j . We also note that our lower bound is 
better than that of Luo [6] for k > 2. 

Proof of Theorem 2.1. For notational simplicity, we presently set Px = fk{x]\), omitting the 
dependence on k and A, and A^; = P^ — Px-i, x = 0, 1, . . .. It is easily seen [H [71 |T0] that the 
probability generating function of Px is 

oo 

g{s) = s'^Px = eM-fc+^+^'+-+«'=), (2.1) 

which implies that 

g'{s) = A(l + 2s + • • • + ks''-^)g{s). (2.2) 

For X > 1, we differentiate (x — 1) times g'{s) and employ the fact that Px = {-^)^-^^ at 
s = to get the recurrence 

k 

xPx = Y,j^P-v ^>1- (2-3) 

We note that (|2.3p is trivially true for x = 0. By definition Px < Pm^ x every x > 0, and 
therefore 

k k 
xPx = Y.^^P-j ^ Y.^^P^K>^ = ^Prn.^Mk + l)/2. 

i=i i=i 

Upon setting X = mfc, A we get m^^A < A/c(A; + l)/2. Therefore m^^A < [Xk{k + l)/2\ since m^, a 
is a non negative integer. 
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As for the left hand side inequahty we note that it is trivially true for k = 1 and A > 0, 
since mi^x = A or A — 1 if A G N, and mi^x = L^J if A N. Therefore we assume that k > 2. 
For < A < 1, the inequality is true since [Xk{k + l)/2j — ''^^^^^ ^ i < g. For A = 1 it is also 
true since e~'^ = Pq = Pi < P2 = 3e~^ /2. Let then A > 1. We will show that Px increases, or, 
equivalently, A^; is positive, for < x < [\k{k + l)/2j — ^^^^^^ + 1. 

From the definition of A,, and (12.11). we obtain 



Ms) = X^5"'A, = (l-sM^). (2.4) 



x=0 



Differentiating h{s) twice we get 

h"is) = A (a ^^s^-' + (^) (A - 2).'^"i + s'^fis)) g{s), (2.5) 

where f{s) = X]j=o«j'S"' is a (A; — 1)*"^ degree polynomial. Next, differentiating x times h"{s) 
and then setting s = 0, we get 

ix+i){x+2) . _ i(i+i) y p 9\p , _L v'^^i n P , 

X — Z^j=l 2 ■^^x+l-j -\ 2 ~ ^l^x-k+1 + Z^j=0 ^J^x-k-j- 

Finally, eliminating successively Px-2k+ii Px-2k: ■ ■ ■ ^Px-k-, by means of (j2.3|) we arrive at 

{x+i)[x+2) ^^^^ ^ (^.^ ^ ^ ^ 1 _ + _ ^ _ 2)P,,_,+i. (2.6) 

Setting = A^-^i-j + Px-j in (|2.6p we obtain 

{x+l){x+2) A ,„ _ /".Y^-L 1^ _L PLllMilli^ A , (A-l)fc(fc+l) 



(2.7) 

Since A > 1, we have Ag = e > and Ai = (A — l)e '^^ > 0. An easy recursion using 
([27711 shows that A^ > for 2 < x + 2 < _^ x ^Iso. This completes the proof of the 

theorem. □ 

Proof of Theorem 2.2. For k = 2, Theorem 2.1 reduces to 3A — 2 < m2^x ^ 3A. Therefore, 
in order to show that m2^x = 3A — 1, it suffices to show that A3A-1 > and A3A < 0. 
However, by (|2.3p . 3A3;j = — 2A3;^_i. Therefore, we will only show A3;^_i > 0. For A = 1, 
A3A_i = A2 = 6-2/2 > 0; for A = 2, A3a_i = A5 = ^ > 0. Let A > 3 and x = 3A - 3. 
Using ()2.6p we have 

(3A-1K3A~2) ^^^_^ ^ _ 3^p^^_^ _ _ 2)p^^_^ ^ _ 3)^^^_^ _ 

By (1131), 

^ ni=i(6A - i)A3A-i = (64A3 - 267A2 + 360A - 156)A3a-7 + 3(A2 + 8A - 12)P3A-8. 

Therefore, A3a_i is positive since A3A-7 > by Theorem 2.1, P3A-8 > by (1.1), and both 
64A^ - 267A2 + 360 A - 156 and A^ + 8A - 12 take positive values. 

For k = 3, Theorem 2.1 reduces to 6A — 5 < 7713 a < 6A. Therefore, in order to show that 
m3^A = 6A — 1, it suffices to show that AgA-j > (1 < j < 4) and AgA < 0. However, 
6AgA = — 5AgA_i — 3AgA_2 bccausc of (|2.3p . We will show then only that AgA_4 > (the 
other three can be treated similarly). For A = 1, AgA_4 = A2 = e~'^/2 > 0. Let A > 2 and 
X = 6A — 6. Using ()2.6p we have 
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(6A-4)(6A-5) 
A 



A6A-4 = (7A - 6)P6A-6 + (8A - 7)P6A-7 " (15A - 12)P6A-8 

= (7A - 6)A6A-6 + (15A - 13)A6A-7 - Pex-s- 



By (ESI) 



nj=4(6A - j)A6A-4 = (1015A3 - 3234A2 + 3396A - 1176)A6a-9 

+(1203A3 - 3610A2 + 3576A - 1176)A6a-io + 2(199A2 - 372A + 168)P6A-ii, 

which is positive, since AgA-g > and Aqa-io > by Theorem 2.1, PgA-ii > by (1.1), and 
their polynomial coefficients take positive values. 

When A: = 4 (/c = 5) we use the same procedure as above to show that Aiqa-j > 
(2 < i < 8) and Aioa-i < (Ai^x-j > (2 < j < 13) and Aisa-i < 0). Therefore, 
m4,A = lOA - 2 (m5,A = 15A - 2). 

Remark 2.1. As k increases the computations become increasingly difficult and lengthy. 
We have used the computer algebra system Derive and a personal computer to check them. 

Remark 2.2. According to the conjecture of Philippou and Saghafi |14) . me, 2 = 39. How- 
ever, by (and ([13])), we presently find that /6(40;2) = 0.0297464817 > 0.0297385179 = 
/6(39; 2). Therefore the conjecture is not true at least for k = 6 and A = 2. □ 



In this note we have derived an upper and a lower bound for the mode(s) of the Poisson 
distribution of order k. Our lower bound is better than that of Luo . We have also established 
the conjecture of Philippou and Saghafi for 2 < k < 5 and A G N, partially solving the 
open problem of Philippou [3[8l[TT]. However, the problem remains open for all other cases. 
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